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Abstract. We study electron transport properties of a monoatomic graphite layer 
(graphene) with different types of disorder at half filling. We show that the trans- 
port properties of the system depend strongly on the symmetry of disorder. We 
find that the localization is ineffective if the randomness preserves one of the chi- 
ral symmetries of the clean Hamiltonian or does not mix valleys. We obtain the 
exact value of minimal conductivity 4e^/7r/i in the case of chiral disorder. For 
long-range disorder (decoupled valleys) , we derive the effective field theory. In the 
case of smooth random potential, it is a symplectic-class sigma model including a 
topological term with 6 — n. As sl consequence, the system is at a quantum crit- 
ical point with a universal value of the conductivity of the order of /h. When 
the effective time reversal symmetry is broken, the symmetry class becomes uni- 
tary, and the conductivity acquires the value characteristic for the quantum Hall 
transition. 

1 Introduction 

Almost 50 years ago, in a pathbreaking paper "Absence of diffusion in certain random lattices" 
[T] Anderson demonstrated that a quantum particle may become localized by a random po- 
tential. In particular, in non-interacting systems of two-dimensional (2D) geometry even weak 
disorder localizes all electronic states [2J31, thus leading to the exactly zero conductivity, a = 0. 

Recent breakthrough in graphene fabrication [4 and subsequent transport experiments 
[5I6I7I8I9] revealed remarkable electronic properties of this material |10lllj . One of the most 
striking experimental observation is the minimal conductivity a of order e^/h observed in un- 
doped samples and staying almost constant in a wide range of temperatures T from 300K down 
to ~ IK. This behavior should be contrasted to well-established results on the conductivity 
of 2D systems where Anderson localization drives the system into insulating state at low T 
[2|3|12|13|14j . An apparently T-independent value of ct ~ e^/h suggests that the system is close 
to a quantum critical point and calls for a theoretical explanation. The aim of this paper is to 
analyze what one should expect for conductivity from the theoretical point of view. We will see 
that, in view of the unconventional character of the graphene spectrum, the nature of disorder 
is crucially important. 

In this paper we consider the conductivity of a single-layer graphene at half filling, e = 0. 
The Drude conductivity [15116] has the value a — Ae'^ /irh at this point. Since this value is of the 
order of conductance quantum, the localization effects should, in principle, become strong at 
half filling. At lowest temperatures T, the system with generic disorder preserving time-reversal 
symmetry falls into the orthogonal symmetry class [12 13 14 implying strong localization regime 
[13j . In this situation, the conductivity at half filling should have a pronounced temperature 
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dependence and get strongly suppressed with lowering T. This conclusion remains valid for 
generic disorder also if the time-reversal symmetry is broken (e.g. by magnetic impurities), so 
that the symmetry is unitary |14| . Remarkably, this is not what is observed in the experiment 
[5|6j . According to what was said above, this might only happen, if at all, for a particular type 
of disorder. 

We identify two broad classes of randomness in graphene — chiral disorder (preserving chiral 
symmetry of clean graphene) and long-range disorder (not mixing the two valleys) — leading to 
the lack of localization and emergence of quantum criticality and associated universal minimal 
conductivity |16ll7j . 



2 The Model 

The tight-binding Hamiltonian of clean graphene is a 4 x 4 matrix operating in the AB space 
of two sublattices (Pauli matrices ai) and in the K-K' space of the valleys (Pauli matrices r^). 
Therefore we introduce the four-component wave function 

^ ^ {4>AK-,4>BK,4>BK' ,4>AK'Y ■ (1) 

In this representation the Hamiltonian has the form 

H = t;oT3crk. (2) 

Here T3 is the third Pauli matrix in the K-K' space and cr = {o'i,(T2} is the two-dimensional 
vector of Pauli matrices in the AB space. The Fermi velocity in graphene is Dq — 10* cm/s. 
In fact, the form of the Hamiltonian ^ is universal and does not rely on the tight-binding 
approximation. The degeneracy of the spectrum in K and K' points is provided by the two- 
dimensional representation of the honeycomb lattice symmetry group, while the expression ([2]) 
is the first-order fc-expansion near these points. 

Let us analyze the symmetries of the clean graphene Hamiltonian ([2]). First, the system is 
obviously uniform and isotropic. Any disorder considered in this paper preserves these symme- 
tries on average, so we do not pay much attention to them here. Second, due to the two valley 
structure of the electron spectrum the whole SU(2) symmetry group exists in an isospin space 
of the valleys. The generators of this group are [12] 

Ax ^ a^Ti, Ay ^ a3T2, A^ ^ aoT^. (3) 

These three operators commute with the Hamiltonian and anticommute with each other. Third, 
there are two more relevant symmetries of the clean Hamiltonian, namely, time inversion oper- 
ation (we denote it To) and chiral symmetry (Co). Using various combinations of Tq, Cq, and 
rotations of the isospin (^o.x.y.z) one can easily construct a whole bunch of twelve symmetry 
operations 

To : A^ aiTiA^aiTi, Cq : A^ -a^ToAasTo, CTq : A ^ - 0-2 n^'^ 0-2X1, 
Tx : Ai-^ 0-2 To ^'^0-2 To, Cx ■■ A^ -ctoTiActoTi, CT^ : A ^ -(TiTqA'^ctiTo, 
Ty : A^ (72T3A^a2T3, Cy : A^ -(JaT2AaQT2, CTy : A ^ -aiT^A^aiT^, 
Tz : A^ (7iT2A^aiT2, Cz ■ A^ -a^T^Aa^r^,, CT^ : A ^ -CT2T2A^CT2T2. 

Note that the transposition changes sign of momentum operators. 

Now we incorporate a disorder of the form VijGiTj in the model, breaking some of the 
above symmetries. The average isotropy of the disordered graphene implies that A^ and Ay 
symmetries of the Hamiltonian are present or absent simultaneously. Below we combine them 
into a single notation A^ , and proceed in the same way with Tl and . In Table [T] we list all 
possible matrix structures of disorder along with their symmetries. 

Generally, the chiral symmetry implies that the Hamiltonian takes block-off-diagonal form 
under a proper unitary transformation. For instance, a generic disorder preserving Cz symmetry 
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Table 1. Symmetries of various disorders in graphene. The first five rows of the table contain disorders 
preserving physical time inversion symmetry Tq. The next four present structures violating Tq. Note 
that each row contains "+" either in the column or in one of the C columns. This means that each 
disorder structure from the first column — when taken alone — produces a finite conductivity at the 
Dirac point. 
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Table 2. Possible types of disorder in graphene leading to a finite minimal conductivity. Each row of the 
table corresponds to a certain universality class of the problem. In the first column we give an example 
of the disorder that belongs to this class (RMF is an external random magnetic field). The relevant 
symmetries are placed in the second column. The symmetry class itself is indicated in the third column 
according to the classification of Altland and Zirnbauer [19120121] . The fourth column displays the 
dominant non-trivial term in the corresponding sigma model. This term is responsible for the absence 
of localization. There are three possibilities [22j: (i) Gade term [24l23j . (ii) Wess-Zumino-Witten term 
[25121] . or (iii) Pruisken 6-term with = tt. The value of the minimal conductivity is given in 
the last column. The first five rows contain chiral-symmetric types of disorder [ripples (o"i,2To) can be 
added to any of them]. They all yield the minimal conductivity 4e^/7r/i (up to power- law corrections in 
a weak disorder strength for chirality, hence "~"), see Sec. |3] The last four rows correspond to the 
case of decoupled valleys (long-range disorder), see Sec. (4] from top to bottom: random Dirac vector 
potential, scalar potential, mass, and any of their combinations. The value of the minimal conductivity 
depends on a particular symmetry of the long-range disorder. 
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can have only off-diagonal matrix elements in the AB space of sublattices [18]. Away from the 
Dirac point (e ^ 0), any chiral symmetry is broken by the diagonal {ctqTq) energy term in the 
action. On the other hand, at half filling the chiral symmetry plays crucial role in establishing 
the transport properties of the system, see Sec. [3l If disorder commutes with Az (long-range), 
the system splits into two copies (decoupled valleys) . The conductivity in this case is discussed 
in Sec. H 

In Table [2] we list all possible situations when the symmetry prevents the localization and 
leads to a finite conductivity at e = 0. This happens if the system possesses chiral symmetry 
or the valleys are decoupled, as we describe in detail below. 
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3 Chiral disorder 

The special class of disorder that we will consider in this section is the randomness that preserves 
one of the chiral symmetries of the clean graphene Hamiltonian |16j . Some possible realizations 
of such type of disorder are listed in Table [2] A peculiar behavior of 2D systems with chiral 
disorder with respect to localization effects has been demonstrated by Gade and Wegner [23|24j . 
They considered a random hopping problem on a square lattice and showed that at zero energy, 
where the system possesses the chiral symmetry, the RG /3-function of the corresponding sigma 
model vanishes to all orders in the inverse conductivity, implying that the conductivity is not 
renormalized. This absence of usual infrared-singular corrections to the conductivity due to 
Cooperon and diffuson loops can be attributed to the fact that the "antilocalizing" interference 
corrections to the density of states cancel the localization corrections to the diffusion coefficient. 
The states at the band center e = are delocalized and the conductivity a{e = 0) takes a finite 
value depending on the disorder strength. According to the classification of Refs. [19|20|21] . 
the system studied in Refs. [23124] belongs to the chiral orthogonal symmetry class BDI if the 
time-reversal symmetry is preserved and to the chiral unitary symmetry class AIII otherwise. 
The peculiarity of the problem we are considering is the Dirac dispersion of carriers. This will 
allow us to prove below a statement that is still stronger than that of Gade and Wegner: we will 
show that for certain types of chiral disorder all disorder-induced contributions to conductivity 
cancel and a{e = 0) takes the universal value Ae'^ /irh. 



3.1 Co chirality 

Let us consider the disorder which preserves the Co-chirality, H = —(i^Hcr^. The random part 
of the Hamiltonian contains matrices cri,2''"3, cri^2Ti^2, and ai^2Ta. According to Ref. [27', the 
random Dirac Hamiltonians preserving the Co chirality and violating the Tq symmetry belong 
to the chiral symmetry class AIII, while the combination of Cq chirality and Tq symmetry drives 
the system into the Bogolyubov - de Gennes symmetry class CI. In both cases, the low-energy 
theory is affected |28I21| by the presence of the Wess-Zumino-Witten term [25] in the action. 

Let us turn to the calculation of the conductivity at half filling, assuming the Cq disorder. 
The conductivity is given by the Kubo formula 



1 

2Tr'h 



■Tr 



(4) 



Here 'Tr' operation implies matrix trace and the spatial integration. Now we use the identity 



0-3 



G«(^)(e)a3 = -G^(«)(-e) 



(5) 



valid when disorder preserves Co-chirality. This allows us to trade all advanced Green functions 
in Eq. (|4]) for retarded ones and thus to present the conductivity at zero energy in terms of 
retarded Green functions. Further, we exploit the following important relation between the 
components of the current operator 



-^^3 = if, 



(6) 



which is the consequence of the Dirac spectrum. At this point, our problem differs from that 
considered by Gade and Wegner [23l24j who dealt with a bipartite square lattice with a non- 
linear electronic spectrum. 

The transformations Eqs. ([5]) and ^ allow us to cast the Kubo formula in the following 
form: 

(7) 



a{e = Q) = -\ 

a—x,y 



At first glance, this expression is zero due to gauge invariance. Indeed, the right-hand side of 
Eq. ([7]) is proportional to the second derivative of the partition function Z[A] = TrlogG^[A] 
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(or, equivalently, first derivative of the current Tr j"G'''^[A]) with respect to the constant vector 
potential A. The gauge invariance imphes that a constant vector potential does not affect 
gauge-invariant quantities like the partition function or the current, so that the derivative is 
zero. This argument is, however, not correct for the zero-order diagram (the one with no disorder 
included) , in view of the quantum anomaly [29l30j . The explicit calculation then yields 



Here S is an infinitesimal imaginary part in the denominator of the Green function. Physically, 
it has a meaning of the electron lifetime or, alternatively, a dephasing rate, and can be thought 
of as modelling processes of escape of electrons in some reservoir or some dephasing mechanism. 
Models with such a uniform constant value of 6 were used in the literature to imitate dephasing 
in quantum dots, see e.g. Ref. [3T]. The corrections to Eq. ([5]) are exponentially small in the 
disorder strength. 

We note that the same universal value of the conductivity in the situation when the only 
type of disorder is the Abelian random vector potential was previously obtained in Ref. [32j ; for 
a certain type of the non- Abelian gauge potential with Gaussian distribution, universal conduc- 
tivity was also obtained in Ref. [33j . The above derivation of the universal conductivity remains 
valid for the case when a magnetic field of an arbitrary strength is applied: the vector potential 
A couples to the current, i.e. to the matrices T3(Ti^2j thus preserving the chiral symmetry. In this 
context, it is worth mentioning the result of Hikami, Shirai, and Wegner [34] who found that 
the longitudinal conductivity in the center of the lowest Landau level of the chiral-disordered 
2D electron gas is equal exactly to cr = Ae'^ /irh in the limit of very strong magnetic field, when 
the Landau level mixing can be neglected. 

3.2 C, chirality 

In this subsection we consider the disorder which preserves the Cz chirality, H = —a^T^Haj^T^. 
The random part of the Hamiltonian may then contain matrices tT3ri^2, cti_2''"3, cti,2''o, and 
CToTi,2- Random Dirac Hamiltonians [27| preserving the Cz chirality and violating the Tq sym- 
metry belong to the chiral unitary symmetry class AIII. WZW terms ^55^ from the two valleys 
cancel out [21j in this case. The combination of Cz chirality and the time reversal invariance 
To corresponds to the chiral orthogonal symmetry class BDI. Finally, the combination of Cz 
chirality and Tz symmetry falls into the chiral symplectic symmetry class GIL In all these 
situations, the resulting theory is of the Gade type [23 24 . 

Let us turn to the conductivity at half filling for a generic disorder preserving the Cz 
chirality. The proof of the universality of the conductivity based on gauge-invariance argument 
does not work now since the Cz chirality transformation of the Green's function generates the 
new vector vertices j'^'^ra instead of physical currents. Nevertheless, for weak disorder we find 
that the conductivity at half filling is still universal, (t(£ = 0) = Ae'^/irh, up to corrections in 
powers of disorder strength (see also [35]). To show this, we first calculate the perturbative 
correction (5ct^^^ to the conductivity of a pure system at the first order in disorder strength and 
find that it vanishes, Sa^^"^ = 0. This implies that the conductivity at the Dirac point does 
not depend on the ultraviolet cut-off of the theory and can be presented as a series in the 
weak disorder strength (indeed the second-order term is found to be finite) [16]. Next, we recall 
that for Cz chirality, the RG /3-function of the Gade- Wegner sigma model |23|24j vanishes to 
all orders, so that there are no singular quantum-interference corrections to a{e = 0) due to 
the soft modes (impurity ladders). This proves that the expansion of (t(£ = 0) in powers of 
disorder converges. Thus for the case of weak disorder the conductivity is universal with small 
corrections in powers of the disorder strength (unlike in the case of the Cq chirality, where the 
corrections are nonperturbative in the disorder strength.) 

Finally, let us mention the case of Cj_ chirality. In an isotropic system considered here, 
both Cx and Cy chiralities are expected to be present simultaneously. This implies that the 
disordered Hamiltonian anticommutes with both ti and T2 and hence is proportional to T3. 




(8) 
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Thus it is split into two equivalent copies. This situation of decoupled valleys corresponds to 
the case of long-range disorder and will be considered in Sec. ID 

3.3 Conductivity at finite frequency or finite temperature 

In this subsection, we analyze the frequency dependence of the conductivity in the case of chiral 
disorder. For completeness, we also keep a small level width 6 introduced above. It was crucial 
for the argument leading to Eq. ^ that the system is exactly at half filling, e = 0. A non-zero 
frequency implies an integration over the energy range of the width which breaks the chiral 
symmetry. When the frequency uj is much smaller than 6, this effect is however negligible, the 
infrared regularization is provided by S, and the universal result ^ survives. In its turn, 6 plays 
no role when w ^ 5: it is the frequency that serves as a dominant infrared cutoff now. Then 
a very interesting new situation arises since to plays a twofold role, leading to two competing 
effects. On one hand, as discussed above, the frequency drives the system away from the chiral- 
symmetric point and thus restores localization. On the other hand, the frequency cuts off the 
singular localization correction. Which of these effects wins? To answer this question, one should 
compare co with the level spacing in the localization area, Z\j(e), where s lu. It turns out 
that A^{s) ^ e ^ UJ for the chiral disorder. This result is rather general and is only based on 
the fact that the operator governing the flow of the system away from criticality couples to the 
energy in the action. Therefore the two competing effects of the frequency (the localization and 
the infrared regularization) "make a draw" - both of them are equally important. The system 
turns out to be, roughly speaking, half way between the chiral fixed point and the conventional 
symmetry. This results in a new universal (frequency-independent) value of the conductivity 

~ /h in the considered regime 5 <^ lu. More precisely, this value depends on the type of 
chirality and the symmetry class of the system away from the Dirac point. 

In the presence of interactions, the temperature T also plays a twofold role, similarly to the 
frequency. On one hand, it induces an averaging over the energy window of the width ~ T, 
thus breaking the chiral symmetry and "switching on" the localization effects. On the other 
hand, the interaction at finite T generates a non-zero dephasing rate t^^{T) cutting off the 
localization corrections. As we showed above, the level spacing A(^{T) is ~ T, so that the result 
of the competition of these two effects depend on the value of Tt^{T). It is natural to assume 
that T^^{T) ^ T [ 13|16j . as follows from the result of Ref. [36] for a two-dimensional diffusive 
metal with Coulomb interaction at cr /h. Then the behavior of the conductivity at low T 
will be qualitatively analogous to that for the case of finite frequency. 

More realistically, one can think about a situation when the disorder is predominantly 
chiral (say, due to ripples and dislocations), but the chiral symmetry is slightly broken, e.g., 
by weak short-range potential disorder. Then the above consideration will be applicable in the 
parametrically broad range of T; at the lowest temperatures, the chirality-breaking effects will 
drive the system into the strong localization regime. 



4 Long-range disorder 

In the previous section we have shown that, if one of the chiral symmetries of clean graphene is 
preserved by disorder, the conductivity at half filling is not affected by localization and is equal 
to Ae'^/'Kh. While various types of randomness in graphene (in particular, dislocations, ripples, 
or strong point-like defects) do belong to the chiral type, the experimentally observed value of 
a is larger by a factor ~ 3, suggesting a different type of criticality. In this section we consider 
another broad class of randomness in graphene — long-range disorder |17j — commuting with 
Az (see Table[T]). We will show the possibility of a new critical state when valleys are decoupled. 
This case has a particular experimental relevance if the conductivity is dominated by charged 
impurities. The ripples [37|38j belong to this type of randomness as well. Numerical simulations 
of graphene with long-range random potential [39l40j provide an evidence in favor of a scale- 
invariant conductivity. 
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The characteristic feature of the long-range disorder is the absence of valley mixing due to 
the lack of scattering with large momentum transfer. This allows us to describe the system in 
terms of a single- valley Dirac Hamiltonian with disorder f32l28!41] , 

H = vo(Tk + a^Vf,{r). (9) 

Here disorder includes random scalar (Vq) and vector (Vi^2) potentials and random mass (V3). 
The Hamiltonian © was considered in Ref. |31] as a model for quantum Hall transition. 

The clean single- valley Hamiltonian ([9]) obeys |32l42j the effective time-reversal invariance 
H — a2H"'"a2- This symmetry (Tl) is not the physical time-reversal symmetry (Tq): the latter 
interchanges the nodes and is of no significance in the absence of inter-node scattering. Below 
we refer to the effective time-reversal symmetry in a single node as the TR symmetry. If the 
only disorder is random scalar potential, the TR invariance is not broken and the system falls 
into the symplectic symmetry class AH, see Refs. |32ll3ll4ll2j . The standard realization of such 
symmetry is a system with spin-orbit coupling; in the present context the role of spin is played 
by the sublattice space. 



4.1 Unitary class 

We start with a more generic case of the unitary symmetry (class A). The TR invariance is 
broken as soon as a (either random or non-random) mass or vector potential is included, in 
addition to the scalar potential. To analyze the critical behavior of the system, we derive the 
nonlinear sigma model for the single- node Dirac fermions [17] at finite energy e. This model is 
formulated in terms of the 4x4 supermatrix field Q operating in Bose-Fermi (BF) superspace 
and bearing the retarded-advanced (RA) structure [31j . The sigma- model action reads [17] 



S[Q] = - Str 



-^(VQf 



'xy 



- ]QV^QVyQ 



Si[Q] + iS2[Q]. 



(10) 



Here 'Str' operation includes the matrix supertrace and the spatial integration. The parameters 
ffxx and axy of the action are the longitudinal and Hall conductivities (in units e"^ /h) given by 
the corresponding Kubo expressions for a single node. In the presence of the constant mass m, 
they acquire the values 



1 

2^ 



l + ^^ii^/(e,m) 



27 

f{x,y) = - 

X 



'xy 



X + y 

arctan 1- arctan 
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m 
"2^ 
x-y 
7 



[fie,m) + f{m,e)], (11) 

(12) 



where 7 is the elastic scattering rate. Note that if the global Tq symmetry involving both the 
valleys is preserved, the observable Hall conductivity is absent due to the cancellation between 
the contributions of the two valleys. 

The imaginary part of the action, iS2[Q], is determined by the well-known topological 
invariant Sti{QVxQ^yQ) = 8i7rN[Q] on the sigma-model manifold ^26^. Possible values of this 
functional are integer multiples of Sni. In Eq. (fTU]) . the topological term is equal to i9N[Q], 
with the angle 9 = 2Traxy + tt. In graphene the mass m is absent, so that axy = 0. Thus, the 



topological angle is 6' = tt. The theory (|T0|) is then exactly on the critical line of the quantum Hall 
transition [26143] , in agreement with the arguments of Ref. 1^ . We thus conclude that graphene 
with a generic (TR-breaking) long-range disorder is driven into the quantum Hall critical point 
(i.e. we have a "quantum Hall effect without Landau levels" [B]), with the conductivity 4cr^. 
The factor 4 here accounts for the spin and valley degeneracy. The value a^j is known to be in 
the range a^j ~ 0.5 — 0.6 e'^/h from numerical simulations [45]. A schematic scaling function 
in this case is shown in Fig. [T^. While formally this conclusion holds for any energy e, in 
reality it only works near half filling; for other e the quantum Hall critical point would only 
be reached for unrealistic temperatures and system sizes. It is worth noting that the sigma 
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Fig. 1. Schematic scaling functions for (a) unitary and (b) symp lectio universality class with topological 
term = n. The ordinary case 8 = is shown by dashed lines. 

model Eq. pU]) is rigorously derived only for a^x ^ 1, whereas axx ^ 1 at e = 0. This should 
not, however, affect the universal critical behavior of the theory governed by symmetry and 
topology of the problem. An alternative route from the Dirac anomaly to a topological term 
employs non-Abelian bosonization, see Refs. |21I41| . 

If a uniform transverse magnetic field is applied, the topological angle 9 becomes energy- 
dependent. However, at zero energy, where axy = 0, its value remains unchanged, 9 = 7:. This 
implies the emergence of the half- integer quantum Hall effect, with a plateau transition point 
at e = 0. An alternative explanation [5] of the half- integer quantum Hall effect is based on the 
concept of the Berry phase tt in graphene [42] . This relates the Berry phase and the anomalous 
topological angle 9 = tt. 

In the above, we have assumed a generic TR-breaking long-range disorder. For specific kinds 
of disorder, additional C or CT symmetries may emerge at zero energy. In particular, if the 
only disorder is a random vector potential (ripples), the system belongs to the chiral class AIII, 
with (7 = Ae^/irh, see Sec. |3l If only the random mass is present, the system falls into class D 
at e = 0. In this case the weak disorder is known to be irrelevant |32I28I41] . implying that the 
conductivity is again Ae^ /nh. 

4.2 Symplectic class 

Let us now turn to the case of preserved TR invariance, describing in particular charged im- 
purities. The system belongs then to the symplectic symmetry class AIL The derivation of the 
(T-model starts with the doubling of ^ variables accounting for the TR symmetry |31j . Then Q 
is a 8 X 8 matrix obeying an additional constraint of charge conjugation Q = Q. 

Since the partition function of the symplectic model is real, the imaginary part of the action 
S2 can take one of the two possible values, or tt. The discreteness of 5*2 suggests that it again 
should be proportional to a topological invariant on the sigma-model manifold. A non-trivial 
topology may arise only in the compact (fermion) sector of Q. The corresponding target space 
is A4f — 0(4n)/0(2n) x 0{2n), where n is the number of fermion species. While for the 
conventional (sufficient for the analysis of conductivity) sigma model n = 1, larger values will 
arise if one considers higher-order products of Green functions. The topological invariant takes 
values from the homotopy group |46j 

^2[Xi=^|„=l] =ZXZ, 7r2[A^F|„>2] =^2. (13) 

The homotopy group in the case n = 1 is richer than for n > 2. Nevertheless, S2 may take 
only two non-equivalent values. Hence only a certain Z2 subgroup of the whole Z x Z comes 
into play as expected: the phase diagram of the theory should not depend on n. Possibility of 
the Z2 topological term in the 2D symplectic sigma model was emphasized by Fendley in Refs. 

[1317]. 

To demonstrate the emerging topology explicitly and to calculate the topological invariant, 
let us analyze the case n = 1 in more detail. The generators of the compact sector are Hermitian 
skew-symmetric 4x4 matrices anticommuting with A = p^: pi/i2, P2fJ-o, P2/^i, and ^2/^3- Here 
Pi and Hi are Pauli matrices in RA and TR space respectively. These generators split into two 
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mutually commuting pairs, each generating a 2-sphere ("difFuson" and "Cooperon" sphere). 
Simultaneous inversion of both spheres leaves Q intact. Hence the compact sector of the model 
is the manifold (52 x S2)/'^2- Thus two topological invariants, Ni^2\Q], counting the covering 
of each sphere, emerge in accordance with Eq. (fT3|) . The most general topological term is 
OiNi +92N2. Due to the TR symmetry, the action is invariant under interchanging the difFuson 
and Cooperon spheres, which yields 61 — 62 = where 9 is either or tt. The explicit calculation 
of the n = 1 topological term yields 

iS2[Q] = ^ Str[(^ ± l)fi2UaUf3] , (14) 



where u = T\/T~^ and Q = T~^AT. The sign ambiguity here does not affect any observables. 
Equation can be cast in the form iS2[Q] = iTT{Ni[Q] + N2[Q]), thus yielding = tt. If 
the TR symmetry is broken, the Cooperon modes are frozen and the manifold is reduced to 
a single diffuson sphere (Q, and hence u, commutes with ^,2) with 182 = i'KNi[Q]. In this 
case, A^i[(5] — (ea/j/SiTr) Str(yl^2MaU/3) = (i/167r) Str(QVa;(3Vy(5/i2). In the case n > 1, the 
topological term is iS2[Q] = iiTN[Q] with N[Q] e I2 = {0,1}, cf. Eq. The statement 

= TV follows from the fact that Q configurations of the n = 1 theory which have an odd 
topological charge N[Q] and yield a negative partition function, are inherited by the n > 1 
theory by virtue of natural embedding of the corresponding manifolds. The explicit form Eq. 
(|14p is, however, not applicable in the n > 1 case, since the manifold no longer decomposes into 
two spheres. 

An ordinary symplectic theory with no topological term exhibits metal- insulator transition 
at a*gp « 1.4 e^/h [48]. If the conductivity is smaller than this critical value, the localization 
drives the system into insulating state, while in the metallic phase, cr > ajp, antilocalization 
occurs. Using the analogy with the quantum Hall transition in the unitary class |26I43| , we argue 
that the topological term with 9 = t: suppresses localization effects when the conductivity is 
small, leading to appearance of a new attracting fixed point at ct^* . The position of the metal- 
insulator transition, crj^, is also affected by the topological term. However, we believe that 
its change is negligible: the instanton correction to the scaling function at large conductivity is 
exponentially small [Mj, and the value of the exponential factor e~^^'^ is still extremely small at 
(T = agp. A plausible scaling of the conductivity in the symplectic case with 9 = tt is sketched 
in Fig. [T}). The existence and position of the new critical point can be verified numerically. 
Recent simulations of graphene '39'40' indeed demonstrate the stability of the conductivity in 
the presence of long-range disorder. Of course, in reality there will be always a weak inter- valley 
scattering, which will establish the localization at lowest T, in agreement with [13|14j . However, 
the approximate quantum criticality will hold in a parametrically broad range of T. 

Before concluding, let us briefly discuss this crossover to localization in graphene with 
charged impurities. The dominant diagonal (cro''o) scattering amplitude does not participate in 
the ultra-violet renormalization [32l28ll3ll6j , in contrast to the case of a finite-range intra- valley 
disorder. The reason is the long-range type of correlations (oc in momentum space) that re- 
sulted in the absence of ultra-violet corrections to the linear conductivity in Refs. [49|39|16|50j . 

1/2 

Near half filling, the scattering rate is 7 ~ hvon-^^p, where rii^p is the concentration of Coulomb 
impurities and the interaction parameter is set to unity. The inter- valley processes occur with 
a much smaller rate 7' ~ a^nimp7, where a is the lattice constant. Note that the value of 7' 
is only weakly corrected due to ultra-violet renormalization, unlike in the model of white-noise 
disorder [13]. Assuming as above t^^ T, the crossover to strong localization (orthogonal 
symmetry class AI) takes place at temperatures of the order of 7'. Using realistic parameters 
from Ref . , we find that the symplectic quantum criticality characterized by a scale-invariant 
conductivity dg.* holds down to 10 — 50 mK. Therefore, the model of long-range charged impu- 
rities is capable of describing simultaneously the two most striking experimental observations 
[516] : the linear density dependence of the conductivity [49139] and the existence of a universal 
minimal conductivity at the Dirac point [17] in a broad temperature range. 
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5 Conclusions 

To summarize, we have studied electron transport properties of a disordered graphcnc layer. 
Wc have shown that the nature of disorder is of crucial importance for the behavior of the 
conductivity. At half filling, the conductivity is of the order of e^/h if the randomness preserves 
one of the chiral symmetries of the clean Hamiltonian or does not mix the valleys. 

For the case of chiral disorder, the exact value of the conductivity is still determined by the 
nature of the infrared cutoff, which may depend on the physical setup. We have analyzed in 
detail the situation when this cutoff is provided by the level width 6 or by the frequency oj. 
In the former case the conductivity takes a universal value Ae^/wh, while in the latter case it 
shows a more complex behavior still being of the order of conductance quantum. 

In the c;ase of the long-range (no inter-valley scattering) disorder, grapli(nic also shows 
quantum criticality at half filling. If the effective TR symmetry of the single-valley system is 
preserved (e.g. when Coulomb scatterers are the dominant disorder), the relevant theory is the 
symplectic sigma model with topological angle 9 = n and the minimal conductivity takes the 
universal value 4c7g* . If the TR symmetry is broken (e.g. by effective random magnetic field 
due to ripples or by external magnetic field), the system falls into the universality class of the 
quantum Hall critical point, with another universal value 4(7^ . 
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